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THE CALCULATION OF DOWNWASH BEHIND SUPERSONIC WINGS 
WITH AN APPLICATION TO TRIANGULAR PLAN FORMS 


By Harvard Lomax, Loma Slxjder, and Max. A. Heaslet 


SUMMARY 

A method is developed consistent with the assumptions of 
small perturbation theory which provides a means of determining 
the downwash behind a wing in supersonic flow for a known 
load distribution. The analysis is based upon the use of 
supersonic doublets which are distributed over the plan form 
and wake of the wing in a manner determined from the wing 
loading. 

The equivalence in subsonic and supersonic flow of the down- 
wash at infinity corresponding to a given load distribution is 
proved. In order to introduce the manipulative techniques 
which are subsequently employed, the unswept wing of infinite 
span is treated for supersonic speeds. The principal applica- 
tion in this report , however , is concerned with the downwash 
behind a triangular wing with leading edges swept back of the 
Mach cone from the vertex. Complete solutions are given for 
the chord plane in the extended vortex wake of the wing and for 
the vertical plane of symmetry. An approximate solution is 
also provided for points in the vicinity of the center line of the 
wake. 

INTRODUCTION 

The linearization of the partial differential equation satis- 
fied by the velocity potential for compressible flow yields, for 
subsonic flight speeds, an elliptic-type equation which is re- 
ducible by means of an elementary transformation to the 
basic equation in incompressible flow. As a consequence of 
this result, wing theory in the subsonic realm employs the 
same concepts and types of analyses that belong to classical 
incompressible theory. At supersonic speeds, the differential 
equation for the velocity potential is hyperbolic in type and 
for wing theory is equivalent to the two-dimensional wave 
equation of physics. In spite of the different character of the 
basic differential equation in the two flight regimes, certain 
formal equivalencies can be set up which are intuitively use- 
ful in the solutions of specific problems. In particular, the 
velocity potentials of a three-dimensional source and of a 
doublet each have analogous forms in the two cases. The 
solution of different boundary-value problems encountered 
in wing theory has been discussed in reference 1, and it has 
been shown how suitable distributions of sources and doublets 
may be used to determine the flow potential associated with a 
given lifting or nonlifting wing. 

The calculation of downwash behind a wing, for incom- 
pressible flow, relies almost exclusively on the use of Prandtl’s 
lifting-line theory which is, in turn, developed from the con- 


cept of a single horseshoe vortex. The conventional approach 
to the general downwash problem is to determine, first, the 
induced field of the simple horseshoe vortex by means of the 
Biot-Savart law and, then, from a knowledge of the span-wise 
distribution of loading over the wing, to calculate finally the 
induced field produced by a vortex sheet composed of super- 
imposed vortices of varying span. 

When downwash calculations are to be extended to the 
case of supersonic wings, it appears at first that the use of 
vortex sheets is inadmissible since no practical equivalent 
to the Biot-Savart law exists. It is, in fact, true that the 
horseshoe vortex no longer plays the outstanding role it 
has at low speeds. However, when a more detailed investi- 
gation is made of the underlying analysis, it becomes apparent 
that vortex theory and the Biot-Savart law can be developed 
from the initial use of a constant distribution of doublets 
over a given surface (e. g., see references 2 and 3). These 
doublets produce a discontinuity in the velocity potential 
at the surface, and, for incompressible theory, the curve 
which bounds the surface can be identified with a vortex 
curve possessing circulation. The proof of the Biot-Savart 
law and the introduction of vortex sheets are direct conse- 
quences of these basic ideas. 

Since, as was shown in reference 1, supersonic boundary- 
value problems involving sources, sinks, and doublets can 
be solved in a manner analogous to that used in low-speed 
theory, a method is therefore provided for an attack on the 
downwash problem for supersonic plan forms through the 
use of doublet distributions. By means of this method the 
downwash immediately back of the trailing edge and at an 
infinite distance behind a wing will be derived and shown to 
agree with the previously published results of P. A. Lager- 
strom (reference 4). 

The present report has three principal aims : First, to out- 
line the theoretical approach to the determination of the 
velocity potential of the flow field associated with a supersonic 
lifting surface and the subsequent calculation of the down- 
wash; second, to apply the theory to the case of a triangular 
wing swept back of the Mach cone and to present the results 
of the complete calculations over the chord plane in the 
extended vortex wake of the wing and on the vertical plane 
of symmetry up to about 40 percent of a semispan; and, 
third, to serve as a guide through some of the more difficult 
mathematical manipulations so that the calculations can be 
extended to other plan forms. A simple first approximation 
is also advanced for the downwash variation about the axis 
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of symmetry, and these approximate values of downwash 
are compared with those obtained from the exact calculations. 

In the theoretical portion of the report, the boundary- 
value problem will be introduced and the solutions, obtained 
from Green’s theorem, will be given for low-speed and super- 
sonic flow. In the section of the report devoted to applica- 
tions, the theory will be used to evaluate the potential func- 
tion at an infinite distance downstream from a lifting wing. 
The theory is next applied to the case of the unswept wing 
of infinite span since the mathematical problems involved 
correspond closely to those for the more general case. From 
this application, a general procedure is developed for treat- 
ing wings with supersonic trailing edges. The final applica- 
tion of the report will be devoted to the triangular wing. In 
all of these applications, it will be seen that the analytic 
expressions which have been obtained in supersonic theory 
for the load distributions over certain plan forms afford a 
means whereby the chordwise distribution of pressure may 
be introduced into the analysis, and, therefore, such expedi- 
ents as lifting-line theory are no longer so essential. 

The entire theory is postulated on the assumptions of 
thin-airfoil or small-perturbation theory and, consequently, 
thickness effects and lifting-plate solutions are additive. For 
the results that are given in the plane of the airfoil, the thick- 
ness effect, which is necessarily symmetrical with respect to 
this plane, is zero. 

The material given in the present report is a combination 
of two previously published NACA Technical Notes (refer- 
ences 5 and 6). 
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LIST OF IMPORTANT SYMBOLS 

velocity of sound in the free stream 
span of wing 
root chord of wing 

complete elliptic integral of the second kind with 
modulus k, k 0 , respectively 


(complete elliptic 
with argument t i 


incomplete elliptic integral of the second kind 
with argument t and modulus k 
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f 2 


dt 
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~EJ_ 

V i So 2 

complete elliptic integral of the first kind with 

modulus k [/£= J n ' V(1 -^y(f=pP)] 

incomplete elliptic integral of the first kind with 
argument t and modulus /<: 


dt 


-t 2 )(l-k 2 t 2 )] 

free-stream Mach number 
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free-stream dynamic pressure pAV'j 


static pressure 
Pl-Pu 


r 

r c 

u, v, w 

A u s 

V 0 

w P 


w w 


Wo 

X, y, Z 

x u yi, Z\ 


Xo 


Vo 


Zo 

a 

P 

do 

Mo 

Po 

$ 


A $ s 


V (X— Xi) 2 +(l/— yi) 2 +(z— zQ 2 
V (x— a:,) 2 — P 2 (y— yi) 2 ~P 2 (s— Si) 2 
perturbation velocity components in the direc- 
tion of the x, y, and z axes, respectively 
u u —u x 

free-stream velocity 

z component of velocity induced by doublet dis- 
tribution over plan form 
z component of velocity induced by doublet dis- 
tribution over wake 
-Fo« 

Cartesian coordinates of an arbitrary point 
Cartesian coordinates of source or doublet 
position 
x 


Co 


Py 

Co 


P_z 


Co 


angle of attack 

Vm,, 2 - i 


P tan 4* 

Mach angle ^arc sin jj-J 

density in free stream 
perturbation velocity potential 


semivertex angle of triangular wing 
sign denoting finite part of integral 


SUBSCRIPTS 

u conditions on upper portion of surface 

l conditions on lower portion of surface 

L. E. conditions at leading edge 

T. E. conitions at trailing edge 

VT wake 

P plan form 

s conditions on discontinuity surface (at Zi = 0) 

I, II, III conditions in regions I, II and III on plan form 
(fig. 4) 

A, B,C,D, I? conditions in regions A, B, C, D, or E, in wake 
of triangular wing (figs. 1 and 2) 


THEORY 

BOUNDARY CONDITIONS 

The proposed problem is one of finding the downwash 
behind a flat plate which supports a loading consistent with 
its angle of attack and plan form. It will be assumed 
throughout the analysis that this load distribution is known. 
Such values were given for several plan forms in reference 7 
and further results can be found in the literature on super- 
sonic wings. 

The load distribution over the wing may be obtained 
from a knowledge of the differences in pressures acting on 
the lower and upper surfaces. Moreover, in thin-airfoil 
theory, where boundary conditions are given in the z=0 
plane (i. e., the plane of the wing), a simple relation exists 
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between local-pressure coefficient and the streamwise com- 
ponent of the perturbation velocity. Thus, assuming that 
the free-stream direction coincides with the positive x axis 
(fig. 1 (a)), and denoting by u the x component of the 
perturbation velocity, it follows that 


Ap__ Pi—p u ^ 2 

2 2 V 0 - 


(u u — u,) 


2A u, 
Vo 


( 1 ) 


where the variables are defined in the table of the symbols. 
Furthermore, from the definition of the perturbation velocity 
potential $ 

$ = J udx (2) 


where a is a point in a region at which the potential is zero. 
Combining equations (1) and (2), the jump in potential in 
the plane 2=0 can be determined by integrating the jump in 
the u induced velocity or, what amounts to the same thing, 
the jump in load coefficient. Thus, 

< 3 > 


downstream of the wing form a discontinuity surface for the 
velocity potential, while A$ s is equal to 0 throughout the 
remaining portion of the xy plane. These conditions, 
together with the fact that the vertical induced velocity w is 
a continuous function at 2=0, are sufficient to determine $ 
throughout space. The values of u, v, and w can then be 
found from the corresponding partial derivatives of $ with 
respect to x, y, and z. The attention in the present report is 
centered on w, the -down wash function. 


SOLUTION TO BOUNDARY-VALUE PROBLEM 


In reference 1, the solutions for boundary-value problems 
of the type under consideration were given for both incom- 
pressible and supersonic theory. The basic differential 
equations satisfied by the perturbation velocity potential are, 
for the two cases, respectively, 


and 


d 2 $ 

bx 2 + dy 2 + bz 2 ~' 0 


2 _ c) 2 <f> _ d 2 <f> _ 

^ dx 2 d y 2 c )z 2 ~~ 


(4) 

(5) 


where the integration extends from the leading edge to the 
point x and A$ s represents the jump in <5 in the xy plane. 
Since load coefficiet A ply must be zero off the wing and since 
u is an odd function in z, the value of u must be zero for 
all points off the wing in the xy plane. It follows that A<f> s 
remains constant at a given span station for all values of x 
beyond the trailing-edge position. 


Vo 





(a) Plan form. 

(b) Sections showing distribution of A<£«. 

Figure 1 .— Sketch showing arbitrary lifting surface together with distribution of A<J? e , the 
jump in perturbation velocity potential in the plane of the surface. 


Figure 1 indicates an arbitrary lifting surface in the 2=0 
plane together with the distribution of A$, for given constant 
values of y and x. In both subsonic and supersonic theory, 
the wing together with the semi-infinite strip extending 


Incompressible theory. — For boundary conditions pre- 
scribed in the 2 = 0 plane; the solution of equation (4) is 


*(*,*,*) 


1 /£>A4>\ , A _, 




where 


( 6 ) 


r= A / (x—x l ) 2 + (y y 0 2 + (2 — 2 i ) 2 


and r is the area for which the integrand does not vanish. 
1 


The terms 


4:wr, 


and 


j_ ( r\ 

4tt Vdsi rj s 


arc equal to the velocity 


potential at x,y,z of a unit source and doublet situated at 
the point x lt y i,0. The remaining terms in the integrand, 
which determine the distribution of source and doublet 
strengths, must be found from known boundary conditions. 
If a lifting- surface fixes the boundary conditions, induced 
vertical velocities on the upper and lower faces of the surface 
are equal so that 

d2j d2j 

and 


Hx,y, Z ) J t Ja$ s i) s dx4yr 


(7) 


Equations (6) and (7) are well known in potential theory 
(reference 3, p. 60), but the derivation usually employs the 
assumption that the value of <3? is zero at all points infinitely 
distant from the wing. This assumption cannot, of course, 
be made in aerodynamic applications where the discontinuity 
surface r extends to x= <» ; as in the case of a lifting wing or 
lifting line with trading vortices. These latter problems, 
with which this report is directly concerned, are of such a 
nature, however, that the induced effects at an infinite dis- 
tance are confined to the plane x= . An investigation of 
the derivation of equation (6) reveals that the conditions 
imposed on in general, can be relaxed sufficiently to permit 
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a discontinuity in a strip of finite width along the entire 
extent of the x axis. The mathematical details of the deriva- 
tion will not be given here but a statement of the restrictions 
on $ at infinity is worthwhile. Thus, denoting by d<t >/dn the 
directional derivative of $ taken normal to a prescribed 
surface, the following conditions apply: 

1. The functions $ and c )$/dn are zero at all points having 
radius vectors which make finite (nonzero) angles with the 
positive x axis, the points lying on a spherical surface of 
infinite radius with center at the wing. (This preserves the 
usual potential theory assumptions except over the portion 
of the spherical surface forming the plane x= °° .) 

2. The values <f> and b$/c)x are bounded at all points 
infinitely distant from the lifting surface and at a noninfinite 
distance from the positive x axis. (This condition places 
restrictions on the values of $ and d<t>/bx in the plane % — 03 .) 

Conditions (1) and (2) are satisfied for a lifting surface of 
finite span, and equation (7) is consequently applicable 
directly to the determination of the velocity potential. As 
an application of the equation, suppose a sheet of horseshoe 
vortices is situated as in figure 2 with bound vortices placed 

X 



Figure 2.— Vortex sheet with bound vortices on y axis and distribution of circulation 

equal to 


on the y axis, trailing vortices extending parallel to the 
positive x axis, and has a spanwise distribution of circulation 

A'h symmetrical to the xz plane and defined for 

Then the velocity potential corresponding to this vortex 
sheet is given by the expression 


$(x,y, 




dx i 

[(*— xi) 2 -f (y— y0 2 +2 2 ] 3/2 


(8) 


portion of the a hVi plane bounded by the leading edge of 
the wing, the lines parallel to the x axis and stemming from 
the lateral tips of the wing, and the trace in the Si — 0 plane 
of the Mach forecone with vertex at the point x, y, z. The 
sign | is to be read “finite part of” and was introduced 

by Hadamard (reference 8) as a manipulative technique 
with the property that 


ITT A(x)dx C x ° A(x)~A(x 0 ) ^ 2 A(x 0 ) 

|Ja (X 0 -X) 3 '*-Ja ( X 0 ~XY « ** (X Q -a) 1 '* ^ 

For purposes of calculations, this was modified in reference 
1 to 


C x o A(x)dx _ fi* A(x)dx 
Ja (x 0 — x) 3/2 Ja \x 0 —x) 3/2 


F(a) — C 


OD 


the asterisk indicating that no upper limit is to be substituted 
into the indefinite integral, the latter being determined as 


F{x) + O 

where 

<?=lim 

x-*xa \_-\/X 0 — X J 

Equation (9) is the direct analogue of equation (6). The 
terms ~ Q-^ s and ~ ^ are equal to the velocity poten- 

tial at x, y, z of a unit supersonic source and doublet situated 
at the point Xi, y u 0, while the remaining terms in the inte- 
grand determine the distribution of source and doublet 
strength and are determined by the known boundary condi- 
tions. 

When the potential function associated with a lifting 
surface is to be evaluated, 


ds! — 5si 

and equation (9) reduces to the form 

Hx,y,z) JJa*. (J: *j s dx,d yi (12) 

In application, the region of integration in equations (7) and 
(12) can be divided into areas occupied, respectively, by 
the plan form and the wake region. Thus, for equation (12), 


When A$ s = constant, a single horseshoe vortex results. 

Supersonic theory. — For supersonic boundary-value prob- 
lems associated with plan forms as indicated in figure 1 (a), 
where the known conditions are given in the z — 0 plane, 
the general solution of equation (5) is given in reference 
1 in the form 




(9) 


where 


T c — V (Z — XiY~ & 2 {y — l/i) 2 — |3 2 (z — 2l) 2 

and the subscript s on the parentheses indicates that the 
function is to be evaluated at The region r is that 


$ (x v 2) = _-^ 2 |7T HZ Alvhy/y, _ 

P{x ’ y ’ z) 2 t I J J pla „ form [(x-Xir-Fiy-ytr-FzT 2 

zflTT- M> s dxidyi /loN 

2x| J J^[(x-x l ) a -0 t (y-yd a -W n { } 

Equation (13) presents a formal solution for the calculation 
of velocity potential and, subsequently, downwash for a 
given surface in terms of A^v Since A3> s was related directly 
to load distribution in equation (3), it is apparent that the 
various known solutions to lifting-surface problems are 
directly applicable. The fact that supersonic theory per- 
mits the determination of load distribution in closed analytic 
form for many simple plan forms provides a distinct advan- 
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tage that is lacking in subsonic theory wherein virtually all 
known results are available only in numerical form. Thus, 
theoretical analysis of problems involving supersonic flight 
speeds can be carried further before recourse to numerical 
methods is necessary. 

APPLICATIONS 

VALUE OF POTENTIAL FUNCTION AT x= o° 

It is possible to show, from equations (7) and (12), that 
the potential functions corresponding to a wing with fixed 
load distribution are identical at x=°° for incompressible 
and supersonic flow. Assuming A$ s known, the values of 
( x , y, z) for the two cases are given, respectively, by the 
equations 


*(*«*) = — ff Jx'dyt 

[ ’ y ’ } 4xJ J plan form [(*-:*)? + (y-y^ + zT^ 

z f f i\i\/lx y (hj y 

4ir J Jwake + 


and 

$(x,y,z) 


—& 2 z \ C C Mt> s dx\dy x 

|J J plan form [ (a: — Xj) 2 — P 2 (y— 2/i) 2 — (3 2 s 2 ] 312 

p 2 z C r A $ s dxidyi 

2 x J J^Kx-xty-fiy-yiy-FzT 2 


for fixed values of y and z the integrals over the plan form 
in both equations approach zero as x increases indefinitely. 
Thus, denoting by x T , B . the value of x x at" the trailing edge 
of the wing, the potential functions at x= °o are given by 
the expressions 


and 


z fV2 

^( m ,y,z). =lim l_ 6/2 A<f>,(x r .B.,yi)dyi 

r- dx i 

J X TB- 2 


-xi) 2 +(y— yiY+z 2 ]*' 2 

ri,/2 

$( m ,y,z) =lim -2^-J_ b/2 a $,(x T . E .,y 1 )dy x 

f* dXy 

J X T-E 


■.E. [(*— Zi) 2 — P 2 (y— yiY— P 2 Z 2 ] 3 


These relations can be integrated once to give for the subsonic 
case 


$(°° ,y,z)= lim 


j: 


6/2 

4ir J-b/2 


A $,(x T . E .,yi)dyi 


( (x — Xj) 

( [(y— yi) 2 +z 2 ] V(x— X,) 2 + (y—y i) 2J rz^x T 

_2_ p/ 2 A<t>,(x T .E.,yi)dy x 
2 ttJ- 6/2 {y—yi) 2 +z 2 


(14a) 


Since, however, A4> s is finite, it follows immediately that 


and for the supersonic case 


9(<^,y,z) =lim 

X — >oo 


— 2 C b l 2 

"X I A4 > s (Xt' B , 

J _(,/ 2 



(x—Xi) 

[ (y — y l) 2 + 2 2 ] V (*■ — Xi) 2 — p 2 (y — 2 / 1 ) 2 — j8V 


* 

X T.E. 


z_ r bl2 A&,(x T .B.,yi)d y l 

2tt J _ 6/ 2 ( y—yiY+z 2 


(14b) 


From these equations, it follows that the sidewash and down- 
wash at x— <*> arc invariant with Mach number, provided the 
load distribution is fixed. In fact, their values depend solely 
on the spanwise load distribution, since the terms correspond- 
ing to the chordwise distribution disappeared in the analysis. 
This has been pointed out elsewhere in the literature. It 
should be stressed, however, that the result which has been 
obtained here states that equal span load distributions in the 
two cases yield equal values of the potential function at 
x= oo . This does not imply that a wing at low and super- 
sonic speeds maintains the same potential function at in- 
finity. When the wing is kept fixed, the distribution of w 
on the wing is fixed, but the load distribution is a function of 
speed. 

DOWNWASH ON AND OFF THE WING 

As a further application, the unswept wing of infinite span 
will be treated for supersonic speeds. In this case, as is 
well known, the induced velocities are zero at all points 
downstream of the upper and lower Mach waves stemming 
from the trailing edge. An abrupt jump in vertical velocity 
therefore occurs at the trailing edge of the wing. Con- 
sideration of this jump for the unswept wing furnishes con- 
siderable insight into the nature of the mathematical diffi- 
culties inherent in the calculation of. down wash on and off 
wings of arbitrary plan form. The calculation for the par- 
ticular case will therefore be followed by a more general 
discussion which will be of value in connection with the later 
treatment of the triangular wing. 


Unswept wing of infinite span. — The pressure distribution 
for the wing of infinite aspect ratio is constant. For this 
so-called Ackeret-type loading, Ap/q is equal to 4 a//3, so 
that, when the leading edge lies along the y axis 



(a) Point P ahead of trailing edge. 


(b) Point P behind trailing edge. 


Figure 3.— Areas of integration for infinite span wing. 


where a is angle of attack. In the wake 


A p 

T 




2aVo 

T 


(15b) 


The downwash, or vertical velocity, will first be found when 
the point is between the Mach waves from the leading and 
trailing edges and then when the point is downstream of the 
trailing-edge wave. 
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Since the wing is of infinite width and experiences no varia- 
tion with y, it is possible to consider the problem at y~ 0. 
Thus, from equation (13), for the case when the point under 
consideration is between the Mach waves from the leading 
and trailing edges with its forecone cutting the wing as shown 
in figure 3(a), 

«->0 7T/3 Jo 

n - a ^(x~X i)2-0 2 22 , 

of* iyi 

h [(x-x^-W-pzY't 

where the symmetry of the problem with respect to the y x 
axis has been used. Computing the finite part of the integral 

*=lim? 




[(a;— Xi) 2 — /3 2 z 2 ] V (x — X\) 2 — /3 2 e 2 — p 2 z 2 


5<t> 5 ,, 2z0V o ae 

w= y—=~r hm 

OZ OZ £ _,0 7T 


r 'x-l3^&+z2 


Jo [(x-xi) 2 —fi 2 z 2 ]-j(x~x 1 ) 2 —p 2 e 2 —p 2 z 2 

Equations (16) and (17) can be evaluated directly to give the 
results 

Vo «/ s 


X; r — /3Z 

y S 


w = — Vo 


For a point behind the trailing-edge wave (fig. 3 (b)), the 
two quantities can be determined in a similar manner. Thus, 


4>=lim 2z P Yoae j c o 

€— *0 IT Jo [ (x X\) 2 ■ 


4" V {x-tyi-pzl 


[ (x — Xi) 2 — j3 2 s 2 ] V (x — aq) 2 — /3 2 z 2 — 0 


__2 ^ r 

-/3V * Jo 1/1 Jc, 


X-P Vn 2 +2 2 (faq 

■o [( x - xO 2 -^ 2 -/^ 2 ] 3 ' 2 


and w is given by the partial derivative of i> withyespect to z. 
The term containing the single integral is zero, since the 
integral itself is bounded for all values of e, while the term 
containing the double integral is readily calculable. Thus, 
the values of the velocity potential and the downwash for a 
point behind the trailing-odge wave are given by the relations 

VoaCo 3 \ 

*"TP|( 

( (20) 

• J 

These results are the familiar equations associated with two- 
dimensional supersonic flatplate theory. 

The point of principal interest in this development is the 
jump in the induced vertical velocity w in passing from a point 
just ahead of the trailing-edge wave to a point immediately 
behind the wave. A study of equations (17) and (19) shows 
that this jump is the result of the discontinuity in the contri- 
bution to the downwash of the term containing the single 
integral. Ahead of the trailing-edge wave, this term yielded 
the result that 

W— — V Q a 

whereas behind the wave, the contribution of the term to w 
was zero. 

The method of attack used in the study of the unswept 
wing can be generalized to apply to arbitrary plan forms. A 
discussion of this case follows. 

Arbitrary plan forms.— -As will be shown later, for any plan 
form with supersonic trailing edge, the jump in the value of 
w in the plane of the wing at the trailing edge can be cal- 
culated directly by means of simple momentum methods. 
At this point, however, it is of more interest to consider in a 
general manner the nature of the integrations involved when 
the point x, y, z is either ahead of or behind the trailing-edge 



(a) Point P ahead of trailing edge. (b) Point P behind trailing edge. 

Figure 4. — Areas of integration for arbitrary plan form with supersonic trailing edge. 


wave. Figures 4 (a) and 4 (b) show a plan form with a straight 
trailing edge with areas of integration indicated for the point 
P in each of the two positions. (The straight trailing edge 
is not a necessary restriction and is only introduced for con- 
venience of notation.) The regions of integration are divided 
under the assumption that the first integration on the plan 
form in equation (13) will be made with respect to y, . When 
the point P is ahead of the trailing-edge wave, therefore, the 
contribution of the wake is zero and the integration over the 
plan form is made to conform with regions I and II. When 
the point is behind the trailing-edge wave, three integrals are 
evaluated corresponding to regions I, II, and III. In the 
case of the infinite aspect ratio, unswept wing region I was, 
of course, nonexistent and, in general, no essential difficulty 
in regard to the limits of integration is introduced by this 
region regardless of where P is situated. In region II, how- 
ever, the problem must be treated in more detail. 

Consider first the case when P is ahead of the wave and 
denote by 4>x Ia the contribution of one side of region II to the 
total potential. Then 
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<t>ii a =lim 

€ — >0 


— 0 2 Z Cx-fl^<?+zt 

2 x Jz 


dx i 


IX 


Agifej/Ocfyi 


»+« [ (*■ — Xi) 2 — 0 2 (y — y x ) 2 — 0 2 z 2 ] 3 ' 2 


where 

Yi =v+-p ^{x—xiY—0 2 z 2 


(21a) 


Similarly, when P is behind the wave and the same subscript 
notation is used to refer to one side of region II, the value of 
$nu is 


$ IIa =lim 

£->0 


— f c ° j„ fpi ^i(x l ,y\)dy l 

2 t J x aXl \ J ll +e[(x-x 1 Y-0 2 (y-y 1 ) 2 -0 2 z 2 ] 3 ' 2 


where Yi is as defined above. 


(21b) 


The contributions of the other side of region II to the po- 
tential will not be considered separately as the behavior is 
identical. When P lies ahead of the wave, e appears in the 
limits for integration with respect to both yi and Xj. This 
corresponds to the situation in equations (16) and (17) and, 
as for that problem, the limiting process is carried out after 
the integration is completed. When P is behind the wave, 
it is not necessary to defer the limiting process, since 




-0 2 z 


2i r 


f C °dxr lim P - 

Jx €->o Jtf+ e [(# 


A4>i {x u y l )dy l 


xi) 2 -0 2 (y-y 1 ) 2 -p 2 z 2 ] 3 ' 2 


and if { } represents the integrand, then 


But since 


o^ 4 ru*-rn* 

(ag 1; yi) 


[ (x - *0 2 - P 2 (y —2/0 2 - 0 2 s 2 } 3/2 ~ 


<M 


with supersonic trailing edges; whereas both the downwash 
and loading are continuous across a subsonic trailing edge. 

TRIANGULAR WING 

Consider a triangular wing (fig. 5) with, leading edges swept 


\ Vo 



Figure 5. — Regions A, B, and C for triangular wing in zoj/o plane. 


for c 0 ^ Xi S x 0 (i. e., { } is bounded for all values of x, in the 
interval of the first integration) ; and further, since 


back of the Mach cone from the vertex. The loading over 
the wing is known to be (references 7 and 9) 


lim M f ! ' +e dw 1 = lim Me = 0 

»o JV «->0 


therefore, for P situated behind the trailing-edge wave, the 
contribution of region Ila is given by 


hr 


— 0 2 z 

2i r 


X>LT' 


Ah^xi ,yx)dy x 


[{x-x l Y-0 2 {y~y i y—0^fi2 


( 22 ) 


The significance of this result is that, when the point P at 
x, y, z is behind the Mach wave from a supersonic trail ing 
edge, the limiting process associated with region II need not 
be considered. When P is ahead of the Mach wave, the 
term e must be retained in the analysis and the limiting 
process used. As was previously noted, the general analysis 
developed in this report places no restriction on the orienta 
tion of the trailing edge; however, it should be pointed out 
that region II exists only for the case in which the trailing 
edge is supersonic. Therefore, the jump in downwash, ob- 
tained from the integration over region II is associated only 


Ap_ A6 0 2 ax, ,, ^ 

2 E n 0^d o V-0% 2 K 

where E 0 is the com plete e lliptic integral of the second kind 
with modulus k 0 =-Jl — d 0 2 and d Q —0 tan \p, \p being the semi- 
vertex angle of the triangle. From equation (3) 


.where 


A3> s =Hsld o W—0 2 yi 2 


H= 


2aVo 

E o 0 


(24) 

(25) 


Setting, for convenience, the velocity potential 

at x, y, z is given by equation (13) to be the sum of the two 
expressions 


^ sHp 2 1 r r V OgX] 2 P 2 y 2 dx t dy l 

2tt IJ Jplanforai [(x— Xi) 2 — 0 2 (y— yd} 2 — 0 2 Z 2 ] 3 ' 2 

* — zH 0 2 \C r YdaW— 0 2 y 2 dx I dy~ 

w 2x J J wakB [( x - Xl Y-0 2 (y- yi Y-0 2 z 2 f> 2 


(26) 


(27) • 
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Equation (26) represents the contribution to the velocity 
potential furnished by the doublets distributed over the 
plan form while equation (27) represents the contribution 
furnished by doublets in the wake. The latter equation is 
the mathematical equivalent in supersonic flow of the sub- 
sonic velocity potential of a sheet of horseshoe vortices 
corresponding to an elliptic span load distribution. Equa- 
tions (14a) and (14b) showed that the expression for 4v at 
x= oo is identical to the velocity potential of the subsonic 
vortex sheet. However, in the vicinity of the x=c 0 line, 
the behavior is entirely different. 

In the present report, equations (26) and (27) will be 
applied to the determination of downwash in the xy and 
xz planes. 

SOLUTION IN THE XY PLANE 

Effect of doublets in the plan form. — For the purpose of 
integration, it is convenient to divide the area behind the 
wing into three regions as shown in figure 5. The division 
lines separating these regions are formed by the Mach cone 
traces from the trailing-edge tips. 

The following symbols will be used in the derivation of 
the expressions for downwash in the xy plane induced by 
the distribution of doublets over the plan form of a triangu- 
lar wing swept behind the Mach cone : 

Ei,E 2 ,Ez complete elliptic integrals of the second kind 
with moduli k x , k 2 , and k 3 , respectively 

K u K 2 , Kz complete elliptic integrals of the first kind with 
moduli k u k 2 , and k 3 , respectively 

*■ VKSO 

*> fe:)C"7) 

*• Vfe) 

71 T - ' 

£(m+m') — V 2 S(/~ m) (m+£) (/— 0 

72 /-• — ■ 

A 1 M + 2 £ 


w+?) + y (g 2 ^ m 2 ) (F-/ 2 ) 

m+m' 

ov +a+ yrr-M 2 ) (?w z ) 

m+V 

£(m + mO + (m + £) (m'~£ ) 

p— m'4 2£ 

M + m' 

P(y—yi) 

c o 


The downwash w P ( x , y, 0) may be obtained by considering 
lim in equation (26) . It can be shown that, in this 

z-»0 02 

case, this limiting process corresponds to taking the partial 
derivative of equation (26) with respect to z and then simply 
setting z equal to zero. Thus the expressions for w P in the 
regions A, B, and G are, respectively, 

... m\^,, (V VQd— ■ ri){y— m) j HP j it rou\ 

Pa ~ 2ir J{ 0 ,1 m (?-V s ) 312 “ V ~ 2 irji 0 lld * (28) 


M ' V (/— i i)(y—i* ) 
* " (e-v 2 ) 3 ' 2 


Wp=- 


P f 

Jj/o+Vo Je 


r - k :~ 

r yo+OgCQ 

-f 3 f Xo +e m+ f ,+9u m 

XTT J »o+^ll J(o 

.. l+«0 ^ 

_H/3 r*o , />' V( m'— v)(v— m) 

2tt J ifp+V o Jm (£ 2 

I+*o 

yp+Vo 

f ‘ « f‘ y -TP=W? i) <*’>+ 

J-yp+Vo Jm — v ) 


-Va+^o 

r i+«o 


ff W-g?) (u — m) 



r 

yo+Vo I -yo+«iFp 

2tt 

P°, 7,d€+ 

Jv o+%io 

f ^ Wfr f 1+ °° /sd? 

J-J's+^iFo I. ho 


L i+9 » 

l+»0 J 


The solution of the three integrals 7), / 2 , and / 3 will be dis- 
cussed in Appendix A. 

The expressions for downwash in regions A, B, and G may 
then be expressed as the following single integrals which 
can be handled by standard numerical methods: 

<3i) 

L l+9 0 

yo+^o 

f Wo / *2 rs 2 (M-M' + 20-2^)l 

J f „ r 2 V«-M)«+/)L £*■ J 


r Nw- 

Ek / (€— m) «+/ 
2£ 2 \ ^2 


82 (m — /T + 2{) — 2 ^m) 


II II !■■■! 
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w P = — 


m 


[J^V 

L i+o 0 


2£(m' — m) 


(K l -E l )dS+ 


Vo+Vo 



-l/o+s^o "| 

f 1+8 ° 1 I w-ri (K 

J{ 0 2fV *3 (Ki— 


Effect of doublets in the wake, 1 — The study of the down- 
wash induced on the xy plane by the doublets distributed 
over the wake will also be divided into the three regions 
indicated in figure 5, and the symbols listed as follows will 
be used in the derivations: 


&A 

1 1 

&B 

5/j / 7 b — 

7 b \S B — v'J 

O'C 

5c \7c v J 

E A , E b- E c 

complete elliptic integrals of the second kind 
with moduli k A , k B , and k c , respectively 

f \ - 

J 


E (h k -) 

incomplete elliptic integrals of the second kind 
>- with arguments l/a A , l/a B , and a c /k c , and 
with moduli k A , k B , and k c , respectively 



k a ,k b ,k c 

complete elliptic integrals of the first kind with 
moduli k A , k B , and k c , respectively 

fe, 


F (h k -) 

incomplete elliptic integrals of the first kind 
► with arguments l/a A , l/a B , and a c /k c , and 
moduli k A , k B , and k c , respectively 



k A 

/S A /y A -v'\ 

\ 7a V 5a — v' ) 

k D 

( 5b TJ(V ( 7b— y'\ 
\v — 7b/ \5 b — v’ ) 

kc 

jVc(v'- 5 c\ 
V 5c Vt c-v'J 

LiL 2 

undefined limits of integration 

x 0 

X 

Co 


2/o 


z 0 


7 A 


7b 


7c 


5a 

5 B 

5 c 

V 


§1 

Co 

' §z 
Co 
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v' — v+2£o 
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r-j-V 


<3(y— 2/1) 

Co 
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/ 2/o—% 
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On integrating equation (27) with respect to Xi and using 
the notation just presented. 




Hp £oZoCo 
2ir ‘ p 


r 

Jl 2 


VO'— v) (17— "') 


(. 2 +z 0 2 ) V£o 2 — . 2 — z 0 2 


dy 


(34) 


The limits on the integral as previously noted differ in the 
regions A, B, and C shown in figure 5; however, in each case 
the limits are roots of one of the two radicals in the integrand. 

It is desirable to express equation (34) in a different form 
in order to obtain an expression for downwash in the plane 
of the airfoil. Integrating by parts. 


$«-= 


H/3c 0 
2tt 13 


Jf" VO'-.) (.-V) tan" 1 - 2 T‘~ 

(L ZoV£o— Z oJ / j 2 


CL, v A rV ’- 2 y 
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£0. 

ZoV£o 2 — y 2 — s 0 2 



(35) 


When lim is considered, it can be shown that in all 

*o-0 &Z 

three regions the contribution to the downwash made by 
the doublets in the wake is given by the expression 


w w = 


m 

2tt£, 


! pi 
0j L 2 


v-\-v f — 2vj 


277 V (v — y) (y — v') 


>s ■yl£o 2 —y 2 dy 


The solution of equation (36) in regions A, B, and G will 
be considered separately. 
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eliminates the linear term in the radical of the integrand, 
and equation (36) becomes 


A. 4~ §4^ 

In region A, L\ = v and Iji=v' . The substitution ij— 


w Wi 


I go 2 7 A /<; \ 

(v—y A ){y A -v') {5a 7a) 


v-ia j j _( $ 

r r+/ i r ii / W-y a 2 ) 
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dt 


(37) 


The identities 
and 


£o 2 = 7a5a 
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are useful in the integration. 

The transformation m=— — - ' 

r — 7 a 

expression for downwasli becomes 


< is next introduced, and the 
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do 
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Integrating the second term by parts and applying the fundamental properties of even and odd functions yields 
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(38) 
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The Jacobian transformation « =sn u reduces the integrals in equation (38) to standard elliptic forms (reference 10), and 


f [<*-••> w-^T vr=u- 


(39) 


Region B 

In region B, Zi — £o and L^—v' and equation (36) may 
be written 


&2 + ,?2 ] Vr(r-+«o++][(&-+ (n-/)] 
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reduce equation (40) to 
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ko 2 (v+v') (S B —y B ) 
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When the transformation u=sn u is made, equation (41) is readily integrable (reference 10), and, after algebraic simpli- 
fication, may be written 
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(42) 


Region C 

In region C, L, = £ 0 and L 2 = — ko, and equation (36) is 
written 


Ww c -- 
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(43) 


The deriviation of the expression for w Wo is similar to 
that for w WjL with the exception that in this case the substi- 


tution 
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is made, and equation (43) may be written 
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SOLUTION IN THE XZ PLANE 


. V (kc 2 d c ~) (1 


=s?>- |> E & 


(44) 


Just as in the study of downwash in the xy plane, so also 
in its study in the xz plane it is convenient to consider 
separately the effects of the doublets distributed over the 
plan form and the wake. The subscript notation for w w 
and w P is the same as before and again w is equal to w w -\-w P . 

Effect of doublets in the plan form.- — In the xz plane two 
regions are indicated in figure 6. Region E lies between the 

Zo 



Figure 6. — Regions D and E in the xqzo plane for a triangular wing. 

Mach wedge from the trailing edge and the line of intersec- 
tion of the two cones from the trailing-edge tips. Region D 
connects this region to infinity. Again the limits of integra- 
tion form the basis of the division into the two regions. 

The symbols listed below will be used in the derivation 
of the expressions for downwash in the xz plane, induced by 
the distribution of doublets over the plan form of the wing. 

Ei,E 5 ,E io E 5o complete elliptic integrals of the second kind 
with moduli /c 4 , k 5 , k,, and /c 5(J , respectively 
Ki,K 5 ,K io ,K 5o complete elliptic integrals of the first kind 
with moduli k it k 5 , k iQ , and k 5Q , respectively 




9<jXi 


V (X—Xi) 2 


k r . 


OvCq 

V (x Co) 2 P 2 Z 2 

-y/ Or — Xi Y 2 — P 2 Z- 
0 o Zi 

VQr— Co) 2 — 1 8 2 z 2 
Co 


Region D 

In region D , equation (26) is written 


T , 


zHp 


rr*s. 


0 


Ve« W-Fyi* 


[ (x — xO 2 — p 2 yi 2 — p 2 z 2 f 2 


dyi 


Integrating with respect to y\ leads to 
zHp j c o ki 


4>f 


{Ki—E^dx i 


it Jo 6 0 xi 

Changing the variable of integration gives 


47 


zHp | 

r \ 9 0 [xkt — V e 0 2 x 2 + P 2 Z 2 (&/ — e n 2 ) \ 

f Ki — E t \ 

' - J 

'o V 9 0 V -\-p 2 z 2 (k ^ — 0 O 2 ) 

W-do V 


dk t 


Taking the partial derivative of C P P[) with respect to z gives 
the expression for the downwash as 


W P =- 


( ki 2 P 2 z 2 ( 

(K, 0 -E ia \ 

1 0o 

( 0oV ' 

\ki 2 -e a 2 ) 

■j8o 2 x 2 +p 2 z 2 (k io 2 -e Q 2 ) 


[xk io — ^e 0 2 x 2 +p 2 z 2 (ki 2 — d 0 2 ) ] + 


x 3 k i d 0 2 


Jo ° ([0 o V+^ 2 s 2 (V-«] 3 

(Ki-EJdk^ 


L ) G 4 2 -0o 2 ) 


( 45 ) 


12 


REPORT 957 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Region E 

In region E, a similar derivation yields the solution 


w P 


t Uvv \i k 5l ' j d n 2 / 
e„ [x — -ylx 2 k^e o 2 -|- p 2 z- (i — k 5 *do 2 ) ] 
V* 2 V^ 2 +/3V(l-V0o 2 ) 

X 


1 Ks-E s do 


x 3 k 5 2 e o 2 


k 5 l—k 5 2 8 Q 2 \[x 2 k & 2 6 a 2 J cP 2 z\l--k' 2 6 a 2 )} 

2 


f 1 , 7 V- ^ <?o ( AkA 2 _ 

Jo W-e 0 2 \[9 0 2 x 2 +p 2 z 2 (k 4 2 - 


2/)_2M3'2 1 

1 


2\ 13/2 


ft. 2 )] 


'4 


(46) 


Effect of doublets in the wake. — The limits of integration 
again necessitate the division of the portion of the plane 
behind the trailing-edge wave into two regions, D and E. 
The following, list of symbols will be used in the derivation 
of downwash in the xz plane, induced by the doublets 
distributed over the vortex wake: 


djy 


&E 


y(x-c 0 y-p 2 z 2 

0Z 

y {x—c 0 y—p 2 z 2 

2 — Co 


0qCq 


«V.= - 


zHp 2 


X' 


0 Q C O 

P 




a ' e 

Jp 2 z 2 +e oV 

Ep, E e 

complete elliptic integrals of the second kind 
with moduli k D and k B) respectively 

E(a' D ,k’ D )\ 

E(a' B ,k' E )J 

incomplete elliptic integrals of the second kind 
with arguments a' D and a' E and moduli k' D 
and k' E , respectively 

k d ,k e 

complete elliptic integrals of the first kind 
with moduli k D and k E , respectively 

F(a' D ,k' D )\ 

F{a' e, k' e) J 

incomplete elliptic integrals of the first kind 
with arguments a' D and a' E and moduli k’ D 
and k' E , respectively 

kp, 

0pCo 

V(2 — C 0 ) 2 — P 2 Z 2 

k' D 

Vl —k D 2 

k E 

V (2 Co) 2 p 2 z 2 

doCo 

k' E 

Vl 

sn 2 {ia D ) 

—a D 2 

snr(ia B ) 

—a E 2 

Region D 


In region D y equation (27) becomes 

f x-p^(y x 2+ 2 2 

8- 

I 

J °0 [ (# 

-x l y-p 2 (y, 2 +z 2 )] 3 ' 2 


Integrating with respect to 2 j, and using the definition of 

Q(\C n 

*l()~ 2 r . n 2 — R 2 ni, 2 

dy 1 


* __zH(x—c 0 ) 




V9 0 V— p 2 yd 


) f 0 = 

-« "Jo (y, 2 +2 2 ) V(2-co) 2 -/3 2 y. 2 -/3 2 ’ 2 ay} 

H(x—c 0 ) goV n VI 


pz^(x—c 0 y—p 2 z i j A , e n V 


( 1 +^ ^i~kn 


- dec 


where 


Jy 1 

d(fio 


Applying the Jacobian transformation u=sn u 


H(x—Co)d 0 2 c 0 2 
D tc@z ■ \k(x — Cn ) 2 ~ ~fi 2 z 2 , 




r° du (48 ) 

Jo 1 +a D 2 k D 2 .sn 2 u 


The expression for <Iv D as given in equation (48) is inte- 
grate (reference 10) and becomes 


4 >, 


H(x—Co)do 2 c 0 2 


_____ in- ( Kp Vl J- a D 2 k d 2 \ / J^\ 

?D 7T/?s V (2 — c o) 2 AjjPz 2 r fl \ k 0 2 ^A+a D 2 )\ia D ) 


iF{a' D , k'j 




uId Vl y&D^ko 2 


Vi+ 


CLd 


iE(a' D , k' D ) —iF(a' D , k' D ) ^]J = - 

Hpz{, ,, , Vl -n'yk'y 

^ ( yi— a V 

A(a/ D , ^ / i) ) E D F(a' D , k' D ) — ^C D £ , (a / B ,^ , i) )^j | 

(49) 

Since it may be shown that 

[Kr,f (a' d, k ' D ) E D F(a' d , k' D ) K B E{a' D , k' d)} = — 

& r (jy- 77 \,.r TS (f <x' p 2 k' P 2 )~] 

V^HbV [} KD -^> a °~ Kd a' D J 

the expression for the downwash w W[> obtained by taking the 
partial derivative of <fv D with respect to z may be written 


_Hp(, E d K D (l—a' D 2 )_ yi—a'o 2 tT ^ w ,_, „ 

w ^— { -¥- d +~~^T d Vl -V ?F ? 1d d ’ d) ~ 


E D F(a' D , k' D ) K D E(a' D , k' />)] 


(50) 
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Region E 

In region E a similar derivation gives 


m j 

1 Vl a' E 2 k' E - 

(Ke-E e , 2 j£ \ 


[ a' E 

t k E 2 a K 


Vl — a' e 2 \EeF{o,' e, k' e )‘ — EeF(cl' E , k' E ) KgE(o/ k' E )] j 

(51) 

Conditions at the trailing edge.- — The value of the vertical 
induced velocity immediately ahead of and behind the 
trailing-edge wave must, of course, be determinable directly 
from equations (21) and (22), respectively, by setting 
T=c 0 ±2(3. If, however, the discussion is restricted to the 
2=0 plane, a much simpler method exists for finding the 
downwash at these points. The approach taken here follows 
essentially that given by Lagers trom in reference 4. 


» 

i Mach lines - 


Vo 


Figure 7. — Sketch of velocity vectors of the air before reaching, on, and after leaving super- 
sonic airfoil. 

Let conditions just ahead of the trailing-edge wave be 
denoted by the subscript 2 and conditions just behind the 
wave by the subscript 3. Figure 7 shows a section of a 
given wing in the plane ?/ = constant. The Mach waves at 
the leading and trailing edge make the angle ji 0 =arc sin 
1/A/o with the 2=0 plane, and the wing is presumed to be 
at angle of attack a. Assuming the trailing edge to be 
normal to the free-stream direction, the variation in the x 
component in velocity when passing through the trailing- 
edge wave can be treated as a two-dimensional problem 
with the condition imposed that u 3 =0 in the 2=0 plane. 

It is known that continuity of flow together with balance 
of tangential momentum across the wave lead to the result 
that the component of velocity tangential to the wave is 
continuous. The tangential components of velocity, V t , 
immediately ahead of and behind the wave are given, 
respectively, by the expressions 

( V t ) 2 ={V Q +u 2 ) cos Mo+«>2 sin mo 
(V t ) 3 =Fo cos M0+W3 sin Mo 

Equating these relations, it follows that 

W 3 = W 2 + U 2 cot Ho 



From equation (23) 

/ AjtA _ 4ff 0 2 ac 0 

\ 2 /* tfolSVW— 0V 


from which it follows that 




e 0 2 \ 

Vtfo 2 — 2/0 V 


(52) 


Approximate values of downwash near center line of 
wake. — The values of downwash which were obtained on the 
xy and xz planes of the wake were exact solutions subjected 
to no restrictions other than those originally imposed by the 
use of the linearized equations of flow. On the center line 
equations (31), (33), (39), and (44) reduce to the considerably 
simpler expressions 


( w ^),-o= 


"0 

M f*» -1 Ki-E l 

it Jo + 


dk x 


(53) 


( w 'c)r=o = 


mv p g.-F. 

,r _Jo &i + 0o 


*0-1 



K 3 E 3 
k 3 \\ +^0^3) 



0 ^),. 0 = 



(55) 


where 


\ _ Hp E c (1 — kc)Kc 

K Ww c)y~ 0 _ 


(^l)ir-O 


dpX\ 
X — Xx 


(56) 


(ks) 


x — x x 

i/=0- eW 


(^^l) y = 0 


00 

X Q — 1 


(*c) 


y=o 


Xp—1 

e 0 ' 


Since the above expressions are relatively simple to com- 
pute, an approximate method based on the generalized 
Taylor’s expansion in the vicinity of the line y=z = 0 can 
be formulated which reduces the tediousness of the calcula- 
tions and gives a good indication of the variation in the 
downwash function in a portion of the wake for points near 
the center line. The next higher terms in the expansion 
can be found without difficulty for the region bounded as 
follows: 


(a) b <V< \ 6 

(b) Both y and 2 lie within the Mach cones 

from the trailing-edge tips 
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The problem resolves itself into one of finding the first 
non vanishing coefficients of y 0 and s 0 in the series 

w/w 0 —A a -\-AiZ 0 -{-Biyo-{-A2Zo 2 -\~C2Zoyo'i~B2y(} 2 ~\~ (57) 


order in y 0 and z B , the downwash function w/w 0 is 


w 

Wo 



Ki — Ej 

ki-\-9 0 



Zg _2 

EgOg \Z 


( 60 ) 


where 


w 0 =— : 


Eg 


277/3 


The value of A 0 is known from equations (53) and (55) to be 


-^oV^+Jo h+do dki ) 


Ag = 


From equations (28) and (36) the expression for total down- 
wash may be written 


WIT 


H/3 2 ( x , 


1 f»a+»o 2/0 — V / (® o— l) 2 — i? 2 . 

: o — 1) J»o — 0 o V W-(2/o-i?) 2 av 


1 % ^ rvo+^i^o £) 

2 J^o- 1 Vs'o-'Vni-f) 


2 /o+»o(^o— €) Ve 0 2 (z 0 —£) 2 —(y 0 —y) 


(e-yT n - 


- drj 


The coefficient B x in the expansion is given by =0 ‘ 

Carrying out the differentiation, with proper regard for the 

singularity in the first integral, it follows that B x = 0. 

Similarly, it can be shown that C 2 — 0 , while the coefficient 

D 1 fd 2 w\ , ,, 

B 2 =h ( \ — 2 ) 1S given by the expression 
2 \oy o /ji 0 =o 


5: 


1 /or/ 2— ^ 

■(2 K *~Y-k2 EA 


)- 


V 

Jo + eo) _ 7?1 (i-^iV 2 “ ^ i— *?] dk> (5 8) 

where the variables have previously been defined. 

In order to calculate the variation with z, it is necessary 
/&$ \ 

to evaluate ^1, = ( % 2 ) where 
\d2oV z o = ‘° 

r M f 3 2/i 2/i (x—Co) 

/3 P Jo 


£ 

77/3 


Vo 

0 „ 2/i f 

arc zy/(x-c 0 r-jv-^z 


5, 2 - 


2/3 


r*’/. 


0 


gJd gX 2 / 3 2 y i 2 

[ (a; — ®i) 2 — /3 z //i 2 — jS 2 s 2 ] 3/2 




The double integral contributes nothing to the coefficient, 
and the remaining portion of the expression can be evaluated, 
without integrating by differentiating twice and using 
Cauchy’s integral theorem 




Thus 


■/!, = - 


1 2 
■EA [2 1 


(59) 


The coefficient .Aa will not be evaluated since the first 
higher order term hi 2 has been found. Thus, to the first 


DISCUSSION 

The variable w\wg (i. e., (wp+w w )/w 0 ) represents the total 
downwash behind the wing divided by the induced vertical 
velocity on the wing itself. If e is the downwash angle and a 
the angle of attack of the wing, then w/w 0 =delda. 



(a) 00=0.4. 

(b) 00=0.0. 

(c) 0 o=O. 8 . 


Figure 8.— Variation of the downwash in the xoVo plane downstream from the trailing edge 
for various span stations. 


In figure 8, the downwash in the xy plane is presented for 
various d 0 ’s and spanwise stations and for all values of x 
from the trailing edge to a point where the asymptotic value 
is closely approached. The region covered in the y direction 
extends from the x axis out to about (1/2) 8 0 , where in the 
coordinate system used x B equals xjcg, y 0 equals /3 y/c 0 , and 
6 0 is the semispan of the wing. Figure 8 can be used to 
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assess the accuracy of the approximation, given by equation 
(60), that the value of downwash is independent of y in 
the neighborhood of the x axis. Within about a half span 
from the trailing edge of the wing no general statement can 
be made as to the variation of w/w 0 in the y direction. For 
distances greater than a half span from the trailing edge, 
however, the variation is quite uniform and w/w 0 deviates 

from its value at y = 0 only slightly for — 7) 0 o . 



Figure 9.— Variation of downwash across the span at various stations downstream of wing 
trailing edge for 0o=O.6. 


For 0o=O.6, a more extensive study was made of the 
variation of downwash with y. Figure 9 represents values 
of w/w 0 across the span for several positions behind the 
trailing edge. Immediately behind the trailing edge the 
value of wjw 0 falls and approaches — °o as the wing tip is 
approached. However, at 0.4 of a root chord behind the 
trailing edge (^ 0 =l-4), wjw 0 rises and reaches the value of 
0.7 as the wing tip is reached. At a; 0 =2.2, the spanwise 
variation of w/w 0 is essentially constant. Although equa- 
tion (60) is applicable only for region A, it is seen from figure 
9 that the approximation that the downwash does not vary 
with y is useful out to about a third of a semispan for all 
values of x. 

The variation of downwash in the xz plane is presented in 
figure 10. The curves represent values of w/w 0 from the 
trailing-edge wave downstream to a point where the asymp- 
totic value is closely approached. In the immediate vicinity 
of the Mach cones from the trailing-edge tips (i. e., r 0 « 1 + 0 O ) 
the curves were not continued because wjw a becomes very 
large and approaches negative infinity as the Mach cone is 
approached. Since this effect results from infinitely large 
values of the radial component of induced velocity at the 
Mach cone, it does not exist in the z Q = 0 plane. Such a 
behavior is consistent with the mathematical idealization of 
infinite pressures at the leading edge and of an abrupt fall 
of load at the trailing edge. However, in an actual flow field 
where these phenomena do not exist the flow will experience 
a milder change in passing across the Mach cone. Even 
in the theoretical results presented in this report the growth 
of the vertical induced velocity in the neighborhood of the 
Mach cone is logarithmic, and the interval in which w/wq is 
appreciably distorted from the general trend is very small. 



(a) 0o=O.4. 

(b) 0o=O.6. 

(c) 0o=O.8. 


Figure 10.— Variation of the downwash in the xozo plane downstream from the trailing edge. 


Some further insight into the behavior of w in the vicinity 
of the Mach cone from the trailing-edge tips can be obtained 
by studying a single vortex which extends infinitely far ahead 
of the origin at an oblique angle to the flow and infinitely 
far behind the origin parallel to the flow (fig. 11)-. The half 
of the vortex which extends ahead makes an angle with the 
free-stream direction less than the Mach angle so that the 


X 



Figure 11. — Induced vertical velocity field for bent and unbend supersonic oblique vortex 
making an angle with the free stream less than the Mach angle. 
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component of free-stream velocity normal to it will be sub- 
sonic. Thus, outside of the Mach cone originating at the 
sudden bend in the vortex at the origin, the flow will be ex- 
actly like that of a linearized compressible subsonic vortex 
with a superimposed uniform velocity parallel to the line of 
the vortex. Inside the Mach cone, however, the flow is 
completely changed. Figure 11 gives an indication of the 
change. The term “bent” vortex refers to the vortex along 
the x axis which is turned suddenly at the origin from the 
angle it had maintained from — °° . The term “unbent” 
vortex on the other hand refers to a vortex which maintains 
the same angle from — °° to + °o . The unbent vortex is 
included in figure 11 for comparative purposes. The figure 
shows that on the 3=0 plane (section A A) the down wash is 
finite and continuous in passing through the Mach cone, but 
that above the 3=0 plane (section BB) the value of w be- 
comes infinite as the cone surface is approached from the 
inside. This behavior at the Mach cone may aid in inter- 
preting the discontinuity in the results for the complete 
wing as given in figure 10. 


Exact linearized theory 

Approximate theory 



(a) t>o=0.4. (b) (c) 0o=O.8. 


Figure 12. — Variation of downwash in the xqzo plane at various positions on xo axis. 

Approximate values of downwash in the xz plane computed 
from equation (60) are compared in figure 12 with the results 
from the exact solution. The approximation that to a first 
order the variation of w/w 0 with z 0 /d 0 is linear with a slope 
— l/E 0 is seen to be useful up to about a third of a semispan. 

Values of wjw a were not computed for points off the xz 
and xy planes; however, the methods given in the report are 
general and directly applicable. The results already given 
would indicate that the approximate solution is valid in the 
vicinity of (1/3)0 O about the x axis. This assumption can be 
checked for large distances behind the trailing edge by con- 
sidering the flow field as x 0 approaches . Thus figure 13 
shows a comparison between the exact value of w/w 0 derived 
by means of the linearized equation and the approximate 
method based upon the use of a generalized Taylor’s expan- 
sion. The agreement is seen to be satisfactory out to about 
one-third of a semispan either vertically or horizontally from 
the x axis. 

Throughout the analysis it was obvious that the calcula- 
tion of the downwash due to the doublets on the wake w w 
was much simpler to perform than the calculation of the 
downwash due to the doublets on the plan form w P . For 
example, the formulas for the downwash on the x axis were 


Exact linearized theory 

Approximate theory 

zo/do 



Figure 13.— Downwash at a large distance behind triangular wing. 


given in terms of the elliptic integrals E and K; for w w the 
evaluation of E and K was sufficient but for w P a numerical 
integration involving E and K was necessary. Therefore, 
in calculating the downwash for wings with plan forms other 



Figure 14.— Variation of the part of downwash on x axis induced by doublets in wake with 
distance downstream in chord lengths, xo=x/co. Triangular wing. 



Figure 15.— Variation of the part of downwash on x axis induced by doublets on the plan 
form with distance downstream in chord lengths, Jo=r/co. Triangular wing. 
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than triangular, it is useful to know in what regions the con- 
tribution of w P to w is small. For this purpose a comparison 
of w w with w P along the x axis of the triangular wing is shown 
in figures 14 and 15. Figure 14 gives the value w w /w 0 , figure 
15 the value of w P /w 0l and figure 16 the total down wash 
(w w -\-w P )/w 0 or just w/w Q . An inspection of figure 15 shows 



Figure 16.— Variation of the total downwash on x axis behind a triangular wing swept be- 
hind Mach cone with distance downstream in chord lengths, xa=x/co. 


that the effect of the doublets on the plan form dies out rap- 
idly behind the point z 0 =l-Mo, that is behind the point of 



Figure 17. — Variation of downwash on x axis with Mach number at various positions down, 
stream of trailing edge. ^=45°. 

intersection of the x axis with the Mach cone from the 
trailing-edge tip. 

An indication of the variation of downwash with Mach 
number is given in figure 17. This figure shows values of 
w/wo on the x axis plotted as a function of M n for various 
values of x 0 . The value of the sweepback angle is 45°, and 
the Mach number range covered could be extended to 1.4 
and the leading edge would remain subsonic. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Nov. 9, 194-8. 


APPENDIX A 


EVALUATION OF SPECIAL INTEGRALS 

Integral /i 

Since there are no singularities in I\, the finite part sign 
may be discarded. The linear term in the radical is elimi- 
nated by the transformation i 7 = (ti + SJ)/( 1-K) and the 
integral becomes 

, /( /- Tl) ( 71 -m) 

Jl “V (Si-TilT! 3 Jm^t, [1-(5,/t0 2 < 2 ] 3/2 • 

«1-M' 

(Al) 


The expressions for Si and y L may be combined to give the 
useful identities 

£ 2 =7i5i 

and 

(Ti — m) ($1 — m') — G“' — Ti) (81 — m) 

Noting that the integrand is an even function, equation 
(Al) may be reduced to the canonical form 


r— o / (m'— Ti)(ti— m) / /— tA f 1 Vl— co 2 7 , Ao , 

il_2 V (81 — Ti)Ti 3 U-//Jo (l-W)*' 2 *" (A2) 

by the substitution 


«i-V , 

co=— > I 

M — Ti 


By the introduction of the Jacobian elliptic functions 
(reference 10) in the transformation w=sn u, the integration 
may be completed, and 


T <■> /(/*'— Ti)(Ti — m) /V— tA f ^ ™ , 

7l - 2 V "157 -tO^ " I sAvJJo 


(5i— T i)ti 

1 / 2€(/- m ) 
rv *i 


(Ki-Ed 


(A3) 


e 

where cd u=^cn u/dn u. 

Integral 7 2 

As the first step in reducing / 2 to canonical form the integral 
is written 


j: [*■ 


P-fOi+MO+M/*' P+fO*+M 7 ) +/*/*' 


2f(i7-f) 


2ffo+0 




drj 


In this case the following identities may be obtained 
directly from the definitions of 72 and S 2 : 


(T 2 — m) (| — 82 ) + (T 2 — 0 (m— 82 ) =0 


and 


V (l+A (m'—v) V (£—17) (i?— m) 

(£ + T2) A' — 8 2 ) + (A — 72) (I + 8 2 ) = 0 


(A4) 


The transformations v = - 2 and &>= — — - t are made, 

1 +t 72— P ’ 

and after algebraic simplifications equation (A4) becomes 


. //-a r. ( a-72) (m+q 0/4-0 (t»-m) a— m) (s-zn r p 
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, p+0 ( -I m , - T2 \~| p 2 

L 2£(T2 + £) \ 5 2 — M/_ J -1 (1 — Ic2 2 03 2 ) 3/2 ^1 — a) 2 


dai- 


By applying the fundamental properties of even and odd 
functions, the first two integrals in equation (A5) can readily 
be integrated. The procedure for handling the finite part 
sign over the third integral will be considered in detail. 
Since, 


J-1 (l-«»)V(l-*2V)(l-«*) dca 2 |1 

dw 


dw 
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then by equation (6) 
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K ’i du 


cn 2 u 


.K^ + lim 


um , — 7= 

*>-» (l-7r 2 2 )Vl 

2 (k-i 5?) 


!3 V2(l — w) (1 — Is 2 )] 

-ijl—lc 2 2 w 2 E 2 ~\ r „ 1 ) 

7 2 ™ V'2(l-A (1 -A 2 ) j - 
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The solution, of equation (A5) becomes, after algebraic 
simplification 

t «i ( / n5 2 (M — m / + 2£) — 2nf \ T jr 


2FV e £2 } 


(A 6) 


Integral I 3 

The procedure for integrating I 3 is similar to that pre- 
viously discussed in connection with I u In this case, the 

integral is canonicalized by the substitution w=-^/ 1 1 and 

the solution may be written 

/3= sW 2 % 7 ~ (A7) 

REFERENCES 

1. Heaslet, Max. A,, and Lomax, Harvard: The Use of Source-Sink 
and Doublet Distributions Extended to the Solution of Arbitrary 


Boundary Value Problems in Supersonic Flow. NACA Rep. 
900, 1948. 

2. von Mises, Richard, and Friedrichs, Kurt O.: Fluid Dynamics. 

Brown University, Providence, R. I., 1942. 

3. Lamb, Horace: Hydrodynamics. Dover Publications (New 

York), 1945. 

4. Lagerstrom, P. A., and Graham, Martha E.: Downwash and 

Sidewash Induced by Three-Dimensional Lifting Wings in 
Supersonic Flow. Douglas Aircraft Co., Inc., Rep. SM-13007, 
Apr. 1947. 

5. Heaslet, Max. A., and Lomax, Harvard: The Calculation of Down- 

wash Behind Supersonic Wings with an Application to Triangular 
Plan Forms. NACA TN 1620, 1948. 

6. Lomax, Harvard, and Sluder, Loma: Downwash in the Vertical and 

Horizontal Planes of Symmetry Behind a Triangular Wing in 
Supersonic Flow. NACA TN 1803, 1949. 

7. Heaslet, Max. A., Lomax, Harvard, and Jones, Arthur L. : Vol- 

terra’s Solution of the Wave Equation as Applied to Three- 
Dimensional Supersonic Airfoil Problems. NACA Rep. 889, 
1947. 

8. Hadamard, Jacques: Lectures on Cauchy’s Problem in Linear 

Partial Differential Equations. Yale University Press, 1923. 

9. Stewart, H. J. : The Lift of a Delta Wing at Supersonic Speeds. 

Quart. App. Math. vol. IV, no. 3, Oct. 1946, pp. 246-254. 

10. Whittaker, E. T., and Watson, G. N. : A Course of Modern 
Analysis. Cambridge University Press, 4th ed., 1940, ch. 22. 


U. S. GOVERNMENT PRINTING OFFICE: 1980 





'v-v; 




AY 



■ l>- - . 


Y- 


, /-, /. 


/. - - 


y - 

"y..- 


, . • •*••. 


-1c. 


-" ,,-r - Positive ■ directions of axes And angies.(forees and morfents^ar^Aown by arrows 


. - ' Axis 1." " ff r. 

Designation ^3*. 

•4^ 

Sym- 

,-bol 

LongrbudinaL— 

Lateral . : . ' 

x i 

. " Y~ 

Normal ^ 

z - 



Force 

(parallel 
to axis) 
symbol 

- Z~ 


X 

■■. Y 

Z 


y Moment about axis - 


Designation 

y :1~ ” 


"Rolling 

Pitching.!— 

Yawing..,.^. 


Sym- 

bol 




L 
M. 
AT 


- Positive 
'direction 


Y *Z - 

Z +X 

Y X yY 


Angle j 


'Designa- 
tion 


Roll— i... 
Pitch...,— 
Yaw.— 


Sym- 

bol 


.-.0 

£ 


if Linear 
^(compo- 
nent, along 
axis)..— 


u 
— v 
-id- 


~^C~ 


Angular 


V 

a 

r 


' ' V 


! C'J. 


Absolute coefficients of moment -y — 

~Y' C G — ^ ' ~~ : G — -~a 

~ l $bS ■-y m ~< l cS> :.^ n ~qbS- - 

(rolling) ^ (pitching) (yawing) 


r Angle of set ofcontrobsurface (relative to neutral „ 
position), 5. (Indicate surface by proper subscript.) 




i 


A 


D - Diameter 
p -Geometric .pitch f 
4 >ID - ' Pitch ratio: v V _ 

, Inflow, velocity YYf 

f Shpstream velocity ' 5 

T^rust, absolute coefricient C T = 


4. PROPELLER SYMBOLS 




f -• 


pv/D* 


- . Q- Torque, absolute coefficient 


~ - ■ / , ■ 

_ . . . , . , . p 

P Power: absolute coefficient 5 (?£= — 5 ^- 

t.. __ , • : J ■ v . w . pn 3 D r ‘ r • 

C, " . Speed-power coefficient == -~vr ’ 

rj i Efficiency -' ;' . , A . . _ - f.p'v : 

' n Revolutions per second, rps - ' % / 


^ _ Q • $ Effective helix angle =tan -1 ( — ) 

Co_ ^D® '•= ~ "" ^ - ■ ■ 7 


Cv~~- 




~ s r" 


y-*- 


1 hp = 76 . 04.kg-m/s'—&50f t-lb /sec 
1 metric ffprsepower= 0.9863 hp “ 
1 mph— 0.4470 mps ' ,jT 
X mps==2.2369 mph 


, - -5. NUMERICAL RELATIONS— AA 


1 lb=0.4536 kg c <. . . . 

1 kg=2.20'16 lb ' ' . : '' -j 
1 mi= 1,609.35 in =5,280 ft 
1 m= 3.2808 it --f" ' 


(- 


Y 


<r . ■ 






